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ANALYSIS OF A DOUBLE FIN-TUBE FLAT CONDENSER -RADIATOR 


AND COMPARISON WITH A CENTRAL FIN-TUBE RADIATOR 
by Henry C, Haller 
Lewis Research Center 


SUMMARY 

An analytical study of a flat condenser -radiator with a double fin-tube 
geometry (closed sandwich) with variable tube side-wall thickness was performed 
for a Rankine space-power electric-generating system. The analysis of the 
double fin radiator included consideration of tube and header pressure drops , 
meteoroid protection for the tubes and headers , along with a detailed presenta- 
tion of the heat rejection analysis and total weight characteristics. The 
double fin-tube radiator is compared to a conventional central fin-tube con- 
figuration on a heat rejection to weight basis for a four -panel radiator con- 
figuration. 

Both fin and tube geometries are compared on the basis of the same power 
level, working fluid temperature, tube and header pressure drop, radiator mate- 
rial, and meteoroid protection criteria. A beryllium radiator for a 1-megawatt 
system and a columbium alloy radiator for a 500-kilowatt system, both at a ra- 
diating temperature of 1700° R, were chosen for the weight and geometry com- 
parisons. 

The conclusion reached indicates a substantial weight savings can be real- 
ized with the double fin-tube arrangement if the tube side-wall thickness can 
be reduced as a result of a possible meteoroid bumper effect of the enclosing 
fins. Weight reductions compared to the central fin-tube geometry of up to 
32 to 39 percent were shown to be possible for the maximum reduction in side- 
wall thickness in the two examples considered. This result further substanti- 
ated the preliminary conclusions given in an earlier reference that compared 
the double and central fin-tube configurations neglecting the effects of 
headers, pressure drops, tube wall temperature drop, and powerplant thermo- 
dynamic cycle considerations. Thus an incentive is offered to further investi- 
gate the meteoroid bumper screen concept and its application in the double fin- 
tube geometry radiator. 


INTRODUCTION 

The generation of large amounts of electric power in space using the 
Rankine cycle concept requires that a large amount of waste heat be rejected 



from the working fluid. This waste heat is the amount of energy that must he 
rejected from the working fluid vapor leaving the turbine in order to complete- 
ly condense it. Since radiation is the only mode of heat transfer for reject- 
ing this energy and since maximum temperatures are limited, the resulting ra- 
diator surface areas and weights are generally large. Previous studies on 
space radiators (refs. 1 to 5) have indicated that a central fin-tube arrange- 
ment is feasible since it reduces radiator weight by reducing the prime surface 
area vulnerable to critical damage by impacting meteoroids. Reference 6 pre- 
sents the analysis and results of a central fin-tube radiator using flat plate 
fins of constant cross section. 

A preliminary comparison of several fin-tube conf igurations was carried 
out in reference 7 that did not include cycle considerations, vapor and liquid 
headers, and tube and header pressure drops. The results of reference 7 indi- 
cated that a substantial weight savings could be obtained by using a double fin 
tube with reduced rectangular tube side -wall thickness instead of the central 
or open sandwich fin-tube arrangements. A reduced tube side-wall thickness can 
result from the assumption that the two fins of the double fin-tube configura- 
tion will act as a bumper screen to incident meteoroids, thus allowing some 
reduction in the thickness of the tube side wall. The double fin-tube configu- 
ration is attractive from a structural viewpoint since it provides a rigid 
struct ure that has a continuous smooth surface that could be the vehicle skin. 
These attractive features motivated the need for a more sophisticated analysis. 

The purpose of this study is to analyze in a more comprehensive manner the 
heat rejection and weight characteristics of a double fin-tube conf iguration 
with variable rectangular tube side-wall thickness, and to more accurately 
identify the potential weight savings over radiators with a central fin-tube 
geometry. The vapor and liquid headers are also included in the geometry ana- 
lyzed since they directly affect the meteoroid protection requirements of the 
radiator tubes and can be a significant portion of the total radiator weight. 
There is also a sizable amount of heat that can be rejected from the vapor 
header for large power systems (ref. l). Working fluid pressure drops were 
also considered for the tubes and headers. 

A rectangular cross-section fin was chosen for the analysis and compari- 
son. A one -dimensional approach was taken in the development of the fin energy 
balance equation with the assumption that the base temperature of the fin is 
equal to the outer surface temperature of the tube armor. The radiator heat 
rejection and weight analysis was carried out for two typical Rankine power- 
planets: a 1-megawatt potassium cycle with a 1700° R beryllium armor and fin 
radiator, and a 500-kilowatt potassium cycle with a 1700° R columbium alloy 
radiator. This report presents the results of the heat rejection and weight 
analysis for the double fin-tube radiator along with a comparison of these re- 
sults with those given in reference 6 for the central fin-tube configuration. 


SYMBOLS 


A surface area, sq ft 
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Ap radiator plant orm area ; sq ft 
vainer able area ; sq ft 


a penetration correction factor 

B constant 

Cp specific heat, Btu/(lh) (°F) 

c a sonic velocity in armor material, y r E a g/p a , ft/sec 
D tube diameter, ft 

E a Young’s modulus of armor material, lh/sq ft 

E angle factor, fraction of the energy leaving a surface that is incident 

upon another surface 

Fyjj vapor header occlusion factor 

g units conversion factor, 32.17 ft/sec 2 

H incident energy, Btu/(hr)(sq ft) 

h heat of condensation, Btu/lb 

J mechanical equivalent of heat, 778 (ft)(lh)/Btu 

K factor, (q-^ + q f ) /eraD^T^ 

Kjj fluid turning loss factor from header to tubes 
k thermal conductivity, Btu/(hr ) (ft ) (°R) 

L minimum half length of fin, equal to L* - Rq, ft 

L* one -half the tube center to center distance, ft 

l actual half length of fin, ft 

N number of radiator tubes 

N c conductance parameter, aZ 2 T-^/kt 

P cycle fluid pressure, Ib/sq ft 

P e powerplant output, kw 

P(0) probability of zero punctures 

Q heat rejection rate, Btu/hr 
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tube radiant heat rejection rate, Btu/hr 

fin radiant heat rejection rate for a fin length 21 
both sides, Btu/hr 

total radiator heat rejection rate, Btu/hr 

vapor header entrance quality 

vapor header heat rejection rate, Btu/hr 

heat rejection per unit length of tube, Btu/hr ft 

radius , ft 

tube side wall to tube centerline dimension, R^ = H Q 
fraction of flow area occupied by one phase 
Reynolds number 
temperature , °R 

vapor saturation temperature at tube inlet, °R 

thickness of fin, ft 

velocity of vapor, ft /sec 

velocity of liquid, ft/sec 

average meteoroid velocity, ft /sec 

weight, lb 

weight flow per tube, lb/sec 
panel width, ft 

normalized distance coordinate, x/z 

fraction of total heat rejected by the vapor header 

fraction of total heat rejected by the tubes and fins 

coordinate measuring distance along lower fin, ft 

normalized distance coordinate, y/Z 

coordinate measuring distance along upper fin, ft 


radiating from 






ft 
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Z radiator tube length, ft 

a,p constants in penetration formula 

5 tube wall thickness, ft 

e surface hemispherical emissivity 

r]* thermal effectiveness 

6 normalized temperature, T/T^ 

[i viscosity, lb/ft-sec 

p density, lb/cu ft (unless otherwise specified) 

a Stefan-Boltzmann constant, 0.173x10“^ Btu/(hr)(sq ft)(°R 4 ) 

t mission exposure time, days 

0,X two-phase -flow parameters 

Subscripts : 
a armor 

b tube base surface 

c tube liner 

cond conduction 
F friction 

f fin 

g vapor phase 

i inside 

LH liquid header 

£ liquid 

m momentum 

o outside 

P particle 

R total 
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rad radiation 

s side -wall 

t tube 

tot total flow, liquid and vapor 
VH vapor header 

X normalized distance coordinate, x/Z 

x coordinate measuring distance along lower fin, ft 
Y normalized distance coordinate, y/Z 

y coordinate measuring distance along upper fin, ft 

0 conditions at tube inlet 

1 base surface 1 

2 base surface 2 

3 radiator inlet conditions 


RADIATOR CONFIGURATION AND THERMODYNAMIC CYCLE 

The general radiator panel configuration considered for the analysis is 
shown in figure 1 where the vapor from the turbine exhaust is distributed to 
the finned tubes by a vapor header. The heat radiated from the vapor header 
and finned tubes causes the vapor to condense. The condensate is then sub- 
cooled and collected in the liquid header before being sent to the condensate 
pump. This scheme can be modified by dividing the radiator into a number of 
nonredundant segments, each of which could be treated as a separate entity. 

The detailed cross-section drawing of the double fin-tube with variable 
tube side wall composed of tube liner inserted in an armor block, which pro- 
vides meteoroid protection, and two rectangular fins is shown in figure 2. The 
tube liner, which is exposed to the working fluid, must be capable of with- 
standing possible corrosion. The liner thickness must also be compatible with 
current fabrication capabilities and structural requirements, and the presence 
of the liner may also substantially reduce the required armor thickness 
(ref. 8). The liner thickness is arbitrarily taken as & c = 0.040 with 
a minimum wall thickness set at 0.020 inch. The liner thickness was increased 
as the inside diameter was increased, so as to provide necessary stiffening and 
strengthening of the radiator tubes. 

For the double fin geometry of figure 2, the liner can be damaged by im- 
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pacting meteoroids in two general ways- The first is by any primary impacts 
occurring on the outer exposed surfaces of the tube block. These impacts are 
assumed to obey the conventional armor penetration and damage relations devel- 
oped for tubes (refs. 8 to 10), with vulnerable area given by 4R^ZN. Accord- 
ingly, the armor thickness 5 a "was determined using the criterion of refer- 
ence 9. The armor thickness, which is a result of the optimization program 
(ref. 6); is applied in full on the upper and lower surface of the tube. A 
second damage source can arise from a spray of particles on the armor block 
side surface ((D 0 - 2t)ZN) resulting from impacts on the fin surfaces. In view 
of the bumper action involved and the obliquity of the secondary impacts, how- 
ever^ a reduction will undoubtedly be allowed in the armor thickness required 
by the tube block side wall to resist the effects of these secondary impacts. 
Since no specific relations are at present available for the determination of 
this side-wall thickness , a parametric variation of the 5 s /s a use ^- to ex- 
amine the effects of reduced side-wall thickness on radiator weight and geom- 
etry. The parameter § s /$ a is defined as the fraction of the armor thickness 
retained on the enclosed side of the tube block. 

The vapor header takes the form of a hollow paraboloid whose wall consists 
of 0.12-inch-thick liner , which was arbitrarily chosen, and whose meteoroid 
armor protection thickness is the same as that required by the tubes. The par- 
abolic shape insures constant velocity in the header. For simplicity, the 
liquid header was designed with a constant diameter and a fluid velocity of 
4 feet per second, so that a very low pressure drop would result. The liner 
for the liquid header follows the same schedule with the inside diameter as do 
the tube liners. However, a maximum liner thickness is set at 0.12 inch. The 
liquid header also has meteoroid armor. 

The assumptions given herein for the double fin-tube radiator geometry, 
with the exception of the fin-tube configuration, also hold for the central fin 
tube that was analyzed in reference 7 and shown in figure 3. 

The thermodynamic cycle used in this analysis is the Rankine cycle, which 
uses a working fluid that undergoes a change of phase. The working fluid is 
condensed in the radiator, which results in a near isothermal condition pre- 
vailing in the tubes and vapor header. In order to show sample results and 
compare the two fin-tube geometries on a heat rejection per unit weight and on 
a radiator geometry basis, two power levels were chosen. Potassium was chosen 
as the working fluid in the cycle for both power levels with a peak turbine 
inlet temperature of 2460° R and a radiator temperature of 1700° R. It was 
also specified that the radiator tubes would subcool the working fluid 100° R. 
Additional cycle requirements such as turbine and generator efficiencies were 
set at 0.75 and 0.90, respectively, with 10 percent of the generator output re- 
quired for accessories and controls. The emittance of the radiator was taken 
to be 0.90, and the effective sink temperature for the radiator was assumed to 
be 0° R. The foregoing values plus the cycle temperatures, working fluid, and 
the power level chosen for the analysis and comparison enabled the determina- 
tion of the total heat rejection rate and the mass flow rate of the working 
fluid for the radiator design inputs. Additional information required from the 
cycle analysis is the quality of the working fluid entering the vapor header. 
The analysis of the thermodynamic cycle used is given in detail in reference 6. 
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HEAT -TRANSFER ANALYSIS 


Approach and Assumptions 

The analysis considers the general case of two rectangular profile fins of 
length l attached to a tube enclosed in an armor block forming a double fin 
configuration as shown in figure 2. Energy input to the fin is comprised of 
heat conduction along the fin from the two tube side -wall surfaces. For any 
specific choice of L/R q the fin length 7 as shown in figure 2 will depend 
on the value of tube side-wall thickness. 

The specific assumptions used in the development of the heat-transfer re- 
lations for the double fin-tube geometry as well as the central fin-tube geome- 
try are 

(1) Radiator emittance is constant with temperature. 

(2) For the determination of radiator temperature variation the radiator 
surfaces act as blackbodies with incident and emitted radiation governed by 
Lambert ! s cosine law. 

(3) Hemispherical radiation to space from both outer surfaces of the radi- 
ator to a 0° R space sink temperature. 

(4) The base surface temperature of the fin is assumed constant along the 
length of the tube and equal to the tube block outer surface temperature. 

(5) Steady-state one -dimensional heat flow exists in the fins and tube 
block. 

(6) Material properties are constant along the length of fin and tube 
block and are evaluated at the fin base temperature. 

(7) The development of the fin and tube angle factors for the radiant in- 
terchange between fin and tube side walls is based on an infinite length of 
tube and fin, with fin thickness assumed negligible in the development. 

(8) The inside tube wall temperature is circumferentially uniform and 
equal to the stagnation temperature of the fluid at the inlet of the header. 


Tube Wall Temperature Drop 


Since the tube block will have thick walls due to the armor required for 
meteoroid protection, a significant temperature drop will occur across the 
block. The assumed heat -transfer paths and the various fin and tube block tem- 
peratures are shown in figure 4. The inside tube wall temperature T* is de- 
termined from reference 6 as 


T* 


1 k h u o \ 

2 Jgh J 


( 1 ) 


8 



where is the radiator fluid stagnation temperature at the header inlet in 

degrees Rankine and Uq is the tube inlet vapor velocity. The turning loss 
factor in equation (l) was set at 1.15 (ref. 6). 

The relation between the tube inside surface temperature and the tube 
block outer surface temperature is based on a simplified approach that assumed 
heat is transferred from the tube inner surface to the exposed surface of the 
tube block by one-dimensional conduction. This heat path was chosen since it 
represents the greatest flow of energy. The heat transmission by conduction is 
assumed to travel a distance 6 a through a cross-sectional area dZ. The 

expression for the heat conduction is 

2kR b 

d %ond =—(**- T b) dZ (2) 

°a 

For simplicity it is assumed that the temperature drop T* - T^ is gov- 
erned primarily by the radiant heat transfer from the exposed surface of the 
block (neglecting the conduction to the fins). The expression for radiation 
from one surface element may be written as 

" 20£R b T b dz < 3 > 

When equations (2) and (3) are combined for d-lcond = dq ra d^ the resultant 
equation for the approximate temperature drop through the tube block wall is 
given as 


ae& a T b + k(T b " T *) = 0 ( 4 ) 

The temperature T-^ obtained from equation (4) is then used in the de- 
velopment of the fin heat -transfer relations. The temperature of the tube side 
wall is also assumed to be equal to the temperature of the tube block exposed 
surface To_ regardless of the thickness of the tube side wall (all values of 


Fin Temperature Profile 


Formulation of equations . - Considering an element of the fin surface in 
figure 2 and employing the previous assumptions y the law of energy conservation 
for an element can be expressed as the energy balance between the net heat 
transfer due to conduction and radiation (ref. 11 ) : 


dQ cond + dQ rad ~ 0 


(5) 


The net internal heat conduction through the element of thickness t and 
length Z is expressed as 


a «cona - s (-»z 


( 6 ) 
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The net radiant heat rejection from "both sides of the element of the fin is 
composed of its emission minus the incident energy. This expression is given 
as 


a «rad = ( 2 < - 


(7) 


The incident energy term H x in equation (7) is composed of the incident en- 
ergy from the tvo side vails to an element of area on the fin 


oT 4 If + F 

°d r dA^-1 dA^-2 


■) 


and the energy leaving an opposing fin surface incident upon the other 

,21 


/ 

J y=o 


cT^ dF 


y ^dA x -dA 


Substitution of expressions (8a) and (8b) into equation (7) yields 

► 27 


^^rad ~ 


- < 1^-1 + ^- 2 ' + 


l 

. -0 


^y ^dAy-dAy 




(8a) 


(8b) 


O) 


Introducing equations (6) and (9) into equation (5) and the dimensionless vari- 
ables 



Y = 


Y 

7 
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JL 


N„ 


oT^Z 2 

"kt 


yield the expression 

% S - ZeHX) - ( Fa A x -l + F dA x -2) - £ (10) 

The angle factors in equation (10) are evaluated using a relation (ref. 

12, eq. (31-58)) that applies to parallel surfaces of infinite length. For the 
configuration of figure 2 the angle factors in question are given as 
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( 11 ) 


dP 






dA x -dA Y 


, ' 2 /B 0 y 

(Y - X)* + 4:1 — ) 


372 


dY 
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*N 

i-i< 

1 - - 

X 


/R 0 \ 2 

9 1 O \ 

X 2 + 4^— j 

l 72 
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( 12 ) 






dA x -2 




2 - X 
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(2 - 


2 / ^° 

x) 2 + 4 t 


172 
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(13) 


Introduction of equations (ll), (12), and (13) and the symmetry of the 
temperature distribution about x = Z, that is, T(x) = T(2Z - x) , into equa- 
tion (10) yield the expression 



where the actual fin length Z is given by the expression (fig. 2) 


Z 


= L + 1 



(15) 
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Computational proced ure . - It was necessary to use numerical techniques to 
solve equation (14) for the temperature profile of the fin of length Z. The 
numerical solutions were carried out on an IBM 7094 electronic digital com- 
puter . 

The second order differential equation d 2 0/d X 2 = B 6^ + f(X) was solved 
for the 6 profile in the interval X = 0 to X = 1. The boundary conditions 
are 0 = 1 at X = 0 and d0/dX = 0 at X = 1, where B is a constant and 
f(X) is a function of X. Kalaba ! s method was used and is now described. 

The term was approximated by a linear function of 6 , namely by the 

first two terms of its Taylor series expansion. Central differences were used 
to obtain an expression for d 2 6>/dX 2 . At each point of a mesh on the interval 
X = 0 to X = 1, the equation d 2 e/dX 2 = B 6^ + f(X) was thus approximated as a 
linear function of 6 . The set of resulting tridiagonal linear equations was 
reduced to two-diagonal form , and these were then solved for the 0 profile by 
backward substitution. An initial G profile guess was used, and each succes- 
sive iteration yielded approximately one decimal place of accuracy. Solutions 
were obtained as a function of the input parameters R c and z/R q . 

' Temperature profile r esults. - Each solution of the fin energy equation, 
which is independent of power and temperature level or tube inside diameter, 
provided a temperature distribution along the fin. Results are plotted in fig- 
ure 5 as a function of position X on the fin for several parametric values of 
N c and z/R q . It is seen that the temperature drop along the fin is very small 
when the conductance parameter E c is small, which indicates a low thermal 
resistance of the fin. Additionally, there is little effect of if R Q when R c 
is small, although a greater temperature drop occurs with increasing if R 0 
ratio. This holds regardless of the choice of the tube wall thickness ratio 
6 s /b a j since only the if R Q ratio is considered in the development of equa- 
tion (14). 


Radiator Effectiveness 


Fin heat rejection . - After the temperature distribution and the slope of 
the temperature distribution curve at X = 0 have been determined, the net 
heat transferred by the fin can be calculated. The rate of heat loss from the 
fin outer surface at any point x on the fin from a differential area NZ dx 
is aeZNT^ dx. The overall rate of heat loss from the pair of fins that 

forms the double fin tube is 


Of = 



eaZNT 4 


dx = 4aZNZT^e 



0 4 (X)dX 


(16) 


This equation is evaluated using the solutions of the fin energy equa- 
tion (14). Comparison of this fin energy rejection to the total heat loss from 
both sides of an isothermal fin-tube section of length 2 (l + R^) can be ex- 
pressed in dimensionless terms as 
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( 17 ) 


where 


* _ 
T) f - 


Qrf 


Hi 


4eaZNZ \ l + — JT- 



e 4 (x)ax 




(18) 


and 0(X) is a function of X for specific values of § s /b a , L/r Q ; and N c . 
Equation (17) can he defined as the thermal effectiveness of the double fin. 


The results obtained from the solutions of equation (17) for fin effec- 
tiveness are presented for an example case in figure 6 as a function of l/R q 
for parametric values of N c and 6 s /& a . In the aforementioned equations de- 
scribing fin effectiveness (eqs. (17) and (18)), the tube wall thickness 5 a 
must be known in order to obtain solutions. In order to obtain the actual 
value of 8 aJJ the tube inside diameter D i; tube liner thickness 5 C , power and 
temperature level; materials; meteoroid protection criteria; tube and header 
pressure drop; and definition of radiator tube vulnerable area must be speci- 
fied. Thus; solutions for equation (17) require complete solutions for the 
entire radiator. Such radiator solutions will be described in detail in the 
section RADIATOR WEIGHT AND GEOMETRY. 


Inspection of the curves in figure 6 reveals that for any fixed value of 
L/R q and S s /b a the fin effectivenss decreases with increasing values of N c - 
The fin effectiveness is also seen to decrease with increasing & s /5 a at con- 
stant N c . When l/r 0 equals zero the value of fin effectiveness does not 

equal zero when the value of Sg/ba is other than 1. This is due to an in- 

crease in the fin length caused by the reduction of tube side-wall thickness 
for any choice of L/R 0 ratio. Thus; at l/R q equal to zero; a fin remains 
of length l = 5 a - & s . 

Tube heat rejection . - The net heat loss from the external surface of the 
tube block is just its radiant emission since there is no incident energy from 

other parts of the system. This energy rejection when compared to the total 

heat loss from both surfaces of an isothermal fin-tube section of length 
2(l + R-fc) can be defined as the thermal effectiveness of the exposed surface 
of the tube. This expression is 


n* 

% 


% 


• Sv 


4oeNZZ {1 + ^)t 4 


1 + 


Rv 


(19) 
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•where z/R^ is defined in terms of R Q , L, 5 a , and 5 s /5 a in equation (18). 
The heat radiation from the tube side -wall surfaces is already included as a 
contribution to the heat radiated from the fin, since in effect the radiant 
interchange acts in a manner similar to the heat conduction down the fin. This 
aspect is' fully explained in reference 11. Solutions for equation (19) also 
require values of armor thickness 5 a and thus are a result of the radiator 
calculation procedure for the double fin-tube block configuration. 

Results of equation (19) for tube effectiveness are shown plotted for an 
example case in figure 7. Tube thermal effectiveness is seen to decrease as 
L/R q is increased, which indicates a decrease in the relative importance of 
the tube portion as the fin length increases. Decreasing the tube wall ratio 
8 g /& a also results in a reduced tube effectiveness, which is caused by a re- 
duction in the tube block outer surface area and thus the energy it can reject 
to space. It is also seen from the figure that for an L/R q ratio equal to 
zero the tube effectiveness does not equal 1 if the ratio & s /& a is less than 1. 
This is brought about because a fin of length 6 a - 5 S remains (see eq. (15)). 

Total fin-tube h eat rejection. - For calculation purposes, it is desirable 
to formulate the total fin-tube thermal effectiveness that can be determined by 
summing the results of equations (17) and (19) for fin effectiveness and tube 
effectiveness, respectively. This expression is 


* Qrf + Ob 1 

^ - / Rb\ 4 ” i 

4aeNZl(l + -fjT-f 1 + ^ 

Figure 8 shows a plot of total fin-tube effectiveness tj* against L/R 0 ratio 
for several choices of conductance parameter N Q and tube block side-wall 
ratio 5 s /6 a for an exaD1 Pl e case. Inspection of the curves shown in figure 8 
reveals that for any fixed value of l/R q and R c the total fin-tube thermal 
effectiveness decreases as the 5 s /b a ratio decreases. This is reasonable 
since as the 5 s /5 a ^atio decreases the amount of isothermal base surface de- 
creases. The fin-tube effectiveness is also reduced by increasing the conduc- 
tance parameter or by increasing the L/R 0 ratio. 

Negligible variation in the magnitude of r}* was observed for the two 
power levels chosen for the comparison. Thus the curves given in figure 8 for 
the 500-kilowatt case also apply for the 1-megawatt power level case. It was 
also noticed that variations of the tube inside diameter had no appreciable ef- 
fect on the value of r]* at a specific value of the L/R q ratio, the 6 s /5 a 
ratio, and the conductance parameter N c . 

PRESSURE DROP CONSIDERATIONS 

Another factor that is required to determine the geometry and weight of a 
radiator is the pressure drop determination in the radiator tubes and headers. 
This aspect of radiator design helps to determine the vapor header geometry and 
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the required tube diameter, the tube length, and the number of tubes. The 
equations presented are for a Rankine cycle condenser -radiator "which assumed 
that vapor flow in the vapor header to be of the same quality as the turbine 
exhaust, two-phase flow in the radiator tubes and all-liquid flow in the liquid 
header. The development of the equations given in this section are given in 
reference 6 . 


Vapor header . - The determination of the pressure drop in the vapor header 
is simplified by assuming that only the gas phase affects the pressure drop. 
This pressure drop is expressed as a ratio of AP to the header inlet pres- 
sure P 3 with the resulting ratio kept constant for comparative purposes. 

This equation, which is for turbulent flow, is 


M _ 0-003S7 p^W 
^3/VH ^3^ e ^*^VH 


( 21 ) 


where Uyjj is the uniform vapor velocity in the parabolic header (based on the 
turbine exhaust quality and neglecting the flow area occupied by the liquid) 
and Re is the vapor Reynolds number based on the vapor header maximum diam- 
eter Pyft- The value of is obtained from the expression 


d vh 


2W Qual^ 
Pgrtu-VH 


, 1/ 2 


( 22 ) 


and the total panel ■width W in equation ( 21 ) is 



(23) 


The term Qualj in equation ( 22 ) is the vapor quality at the entrance to the 
vapor header. The number of condensing tubes R in the previous equation is 
determined from the tube pressure drop analysis in conjunction with the opti- 
mization procedure used for this fin-tube geometry. 


The amount of heat rejected to space from the parabolic vapor header is 
determined from the following expression: 


OvH = ^re«3 X VH = 2 -09 


(24) 


where the factor Fyg is defined as the vapor header occlusion factor for 
radiant emission to space and given a value of 0.85 for this analysis (ref. 

13). The analysis and results of reference 13 (fig. 4) although primarily for 
tubes and fins are assumed valid for the occlusion of the header-panel arrange- 
ment shown in figure 1 . 


Radiator Tubes 

The pressure drop in the radiator tubes where flowing vapor is condensing 
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was computed from a combination of several basic flow and energy equations that 
are given in detail in reference 6. The flow model (ref. 14) used assumed that 
at any given section perpendicular to the flow direction, the temperature and 
pressure in both the liquid and vapor are uniform and the same for both phases. 
This flow model also assumed turbulent flow with liquid and vapor velocities 
uniform in each phase at a given cross section, but that the two velocities 
were not necessarily equal. Pressure drops were computed for a series of in- 
cremental tube lengths, and the pressure drop for the whole tube was obtained 
by summing the incremental drops. 


The total change in pressure for an entire radiator tube is comprised of a 
frictional and a momentum component. The friction pressure drop is described 
by the expression 


dP F 




1.8 

S 


/O- 092 Pg 

\ SD i P g 



(25) 


where Re = ^Wg/nD^g and is a function of X. The differential form of 

the change in pressure due to a change in momentum is 


- dP m = 


1_ 

gA 


+ 4 (»g A g u2 )_ 


(26) 


A third relation is required that relates the increment of tube length and 
the increment of condensate formed: 


KecrjtDjT^ 
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where for the double fin-tube radiator the definition of K is 


K = 




(27) 


(28) 


Equations (25), (26), and (27) along with equation (l) are solved simultaneous- 
ly for dW, dP m , and dPji. The total change in pressure for each increment can 
then be found from the relation 


dP = dP m + dP F 


( 29 ) 


and the total pressure change for the entire radiator tube can be found by 


16 



I 


summing the incremental changes. For comparative purposes, the total tube 
pressure drop is expressed, as the ratio of pressure drop to the tube inlet 
pressure . 

The pressure drop associated -with turning the vapor from the vapor header 
into the radiator tubes and accelerating the flow can be calculated from the 
expression 

1 u 0 

^entrance = % 2 p g T (30) 

where the tube entrance loss factor Ktt is given the same value used in equa- 
tion (l). 


Liquid Header 


The pressure drop in the liquid header is obtained by applying Fanning’s 
equation with a friction factor for turbulent flow. This expression is 


0.00102 PjjV^W 

DL H (R e P °‘ 2 


( 31 ) 


where Re is the Reynolds number corresponding to the maximum liquid velocity 
V LH that, in these calculations, was taken as 4 feet per second. The liquid 
header diameter D LH , which is assumed constant, is determined by applying the 
continuity equation at the header exit to give the expression 


^LH “ 


W 


, 1/2 


C P £ V LH, 


( 32 ) 


RADIATOR WEIGHT AND GEOMETRY 

The heat-rejection analysis of the fin and the tube was nondimens ionali zed 
so that the results could be used for general design purposes. Using the pre- 
vious results for analyzing the merits of the variable tube wall double fin- 
tube configuration in radiator designs is impractical without consideration of 
the total weight of the vapor and liquid headers as well as the tubes and fins. 
It is necessary, therefore, to consider the ratio of heat rejection per unit 

weight O ./¥ and the influencing effects of radiator geometry limitations on 
u? e j 

the maximum heat rejection per unit weight. 


Armor Thickness 

In order to determine both radiator weight and geometry, the effects of 
meteoroid penetration must be considered on both the tubes and headers. This 
will dictate the required armor protection thickness needed for the radiator. 


I 
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The tube armor thickness 6 a is determined using the meteoroid protection cri- 
teria given in reference 9, which is based on a comprehensive appraisal of the 
available data and theories concerning the meteoroid penetration phenomenon. 
According to reference 9, the resultant equation for the armor thickness 6 a 
is given by the expression 



•where in the previous equation 

a =1.75 

(3 = 1.34 

Pp = 0.44 g/cu cm 

Vp = 98,400 ft /sec 

a = 0.53X10" 10 g*V(sq ft) (day) 

Insertion of these constants into equation (33) along with utilizing Young’s 
modulus in the definition of sonic velocity of the material yields the more 
compact form 



The total exposed area to be protected by direct impacts is assumed 

to be the outer surface of the vapor header and the projected area of the tube 
block. The liquid header contribution is assumed to be negligible since its 
surface area is small compared to that of the vapor header. Thus 

Av = At + Ayg (35) 

The radiator tube projected area is given by the expression 


h = = 


Qrejlf 


a£T t^ + r 


where l /R^ is obtained from equation (18) and t}^ from equation (20). The 
fraction of the total radiator heat that is rejected by the tubes and fins 
is defined by the expression 
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The vulnerable area of the vapor header is assumed to be its full surface 
area and is given as 


vH 


Ore.j X VH 

0eF VH T 3 


(38) 


•where Fy^ is the occlusion factor for the effect of the radiator panel on the 
heat rejection of the vapor header surface. 


Combining equations (35) to (38) yields 

1 “ X VH 


a = 
v oe 


u/H 




(39) 


Weight Ratio 

The total heat rejection per unit weight of a fin-tube radiator can be ex- 
pressed as 


^rej 

W 




ej 


NZ 


W 

NZ 


(40) 


where for radiation from both sides of the radiator and from the vapor header 
Qrej 


NZ 


= oe 


^ot^R + R q ) 


2.09 


NZ 


(41) 


Insertion of equation (23) into equation (41) yields the following form of the 
total radiator heat rejection per unit length of tube: 


Or 


NZ 1 = oeR o 


(i * £)(*« * — ^ 


(42) 


where the outside tube diameter (block width) D Q can be expressed in terms of 
the liner thickness 6 C , the inside tube diameter and the armor thickness 
6 a and can be given as 


D = D. +25 + 26 

o i c a 


(43) 


The total weight of the radiator is comprised of the individual weights of 
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the vapor header, the liquid header, and the tube -fin panel. The veight of the 
liquid inventory in the subcooler portion of the radiator tubes is neglected 
since the subcooler length is small. The vapor header -weight is given by 
(ref. 6) 


¥i 


VH 


- T * [( d vh 


+ 5- 


VH„ ) 5 VH c Pc + ( i Vh 


+ 2& VB c + 5 a 


^ 3 a p a 


(44) 


and that of the liquid header and condensate by 




2 

£ 4 LH + p c 6 LH c ( D LH + & LH e ) + Pa 6 a( D LH + S a + 25 LH C 


(45) 


The fin and tube panel weight can be calculated using the expression 


W tf = 4Ztp f NZ + p c | WZ[(D 1 + 2S C ) 2 - D?] + HZp^^ - £ (D ± + 2& c ) 2 ] (46) 


When t = aZ 2 T b /kW c , D 0 from equation (43), 7 from equation (15), and 

R b = + 5 C + (S s /& a )& a are introduced into equation (46), the equation for 

the tube and fin panel weight becomes 
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- Jt(R. + & c ) 2 > (47) 


The total radiator weight is then obtained by summing the results of equa- 
tions (44), (45), and (47). The denominator of equation (40) can be found by 
dividing the total radiator weight by the total tube length NZ. This result 
along with the results of equation (42) when inserted into equation (40) yield 
the radiator heat rejection per unit weight. The peak value of Q^j/W can 
then be obtained by plotting the results of equation (40) as a function of 
L/R for specific values of conductance parameter N c and tube side-wall ra- 
tio 0 s s /s a . 


Panel Geometry 

In addition to the important aspect of minimizing weight for practical ra- 
diator designs, it is also of interest in most cases to investigate the geom- 
etry of the radiator as it might affect the integrating of the vehicle and ra- 
diator. Planform area, aspect ratio, and fin thickness are three facets of the 
geometry of the radiator panel that must be determined in order to satisfy 
radiator -space vehicle integration and structural and fabricational require- 
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ments of the fin-tube configuration. 


Radiator planform area Ap is obtained from the equation 


A p = 2NZ(L + R q ) = 

2 


( 48 ) 


It is seen from this equation that planform area -will vary inversely with over- 
all fin-tube thermal effectiveness rjp for a specific choice of power and tem- 
perature level. The planform area will generally increase with increasing 
L/R q because decreases as the l/R q ratio is increased (fig. 8). 

Another interesting factor with respect to the geometry of the radiator is 
the magnitude of the fin thickness. Radiator applications might require that 
the fin have structural or fabricational qualities that could result in non- 
optimum weights and dimensions. Fin thickness t for the closed sandwich con- 
figuration can be calculated from the expression 


t 


orh 2 

-ST 


(49) 


The panel aspect ratio, which is defined as the ratio of panel width ¥ 
to tube length Z, is obtained by using equation (23) for W and the results 
of the pressure drop calculations for Z. 


Method of Solution 

Simultaneous solution of equations (4), (20), (21), (29), (33), and (39) 
requires inputs of inside tube diameter Dj_, tube and vapor header pressure 
drop ratios, power level, temperature level, and the properties of the mate- 
rials and cycle fluid. This results in values of the parameters 5 a , T-^, Xyjj, 
N, and uq for a selected variation of L/R 0 , conductance parameter N c , and 
the tube side-wall ratio 3 s /6 a . 

Important results required for weight and geometry calculations will in- 
clude the number of tubes N, tube length Z, tube outside radius R Q , panel 
width ¥, and the inside diameter of the vapor header D-^. An electronic 
digital computer was used to obtain the desired results that required an itera- 
tive type solution. 


RESULTS AND DISCUSSION 
Calculation Inputs 

Calculations making use of the resultant equations developed in the analy- 
sis for the double fin-tube geometry require inputs such as inside tube diam- 
eter, radiator vapor inlet temperature, cycle power level and conditions, mate- 
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rials of construction, meteoroid protection criteria, tube block side-wall 
thickness ratio, and pressure drop in the tubes and header in order to com- 
pletely specify a radiator solution. For this reason two specific cases are 
used in order to show the effects of some of the previously mentioned variables 
on heat rejection per unit weight and on radiator geometry. 

The first case considered is a 1-megawatt electrical output powerplant 
with the radiator at 1700° R. The tube armor and the fin were assumed to be 
made of beryllium, and the tube liner was assumed to be a columbium alloy. The 
second case is a 500-kilowatt system with the radiator at 1700° R. For this 
case, the tube armor, tube liner, and the fin were all taken to be columbium - 
1-percent zirconium alloy. These two cases both used a peak cycle temperature 
of 2460° R and potassium as the cycle fluid. Inside tube diameters of 3/8, 
l/2, 5/8, 3/4, and 1 inch were chosen with tube lengths increased to allow 
100° R of subcooling. A 500-day mission time and a probability of no puncture 
P(0) of 0.995 were chosen for the calculation of meteoroid protection thick- 
ness. Pressure drop ratios for the two cases were set at AP/P = 0.02 for the 
vapor header and AP/P = 0.05 for the radiator tubes. The emittance of the 
surface coating on the fins, tubes, and headers was taken to be 0.90. Radiator 
material properties were assumed constant with temperature and evaluated at 
1700° R. The all columbium alloy radiator had a density of 530 pounds per cu- 
bic foot, thermal conductivity of 34 Btu per hour per foot per °F, and a modu- 
lus of elasticity of 0. 202xl0 10 pounds per square foot. The beryllium radiator 
used a density of 115 pounds per cubic foot, a thermal conductivity of 51.5 Btu 
per hour per foot per °F, and a modulus of elasticity of 0.397X10 10 pounds per 
square foot. The beryllium radiator used a columbium liner with the material 
properties previously mentioned. 

Calculations using the given inputs and specifications were also made for 
the central fin-tube geometry that used the method described in reference 6. 


Radiator Weight 

Radiator heat rejection per unit weight Q^j/W was plotted for each 
value of inside tube diameter and tube block side-wall ratio 6 s /b a 

chosen for the comparison over a range of the parameters N c a,nd L/R q . Re- 
sults showing the variation in heat rejection rate per unit weight as a func- 
tion of L/R q for several values of conductance parameter N c are shown 
plotted for two sample cases in figure 9. The figure shows results for the tube 
block side-wall thickness to tube armor thickness ratio 5 s /5 a equal to 0.5 
with tube inside diameters of 3/8 and 5/8 inch chosen for the 500-kilowatt and 
1-megawatt cases, respectively. Each constant N c curve is seen to peak at a 
specific value of L/R q with the value of the l/r o at peak Q re j/ ¥ increas- 
ing as N c is increased. 

Plotting the maxima for each N c curve of figure 9 and the additional re- 
sults for the 5 s /6 a = 0 and 1.0 cases yield a performance map (fig. 10) that 
plots peak heat rejection per unit weight against L/R . It is seen from fig- 
ure 10 that decreasing the armor block side-wall thickness results in a sub- 
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stantial increase in the value of Q^gj/w. According to figure 10, the maxi- 
mum Qrej/ W occurs at an N c approximately equal to 0.75 for the 500-kilowatt 
system at a tube inside diameter of 3/8 inch. The results for the 1-megawatt 
system indicate this maximum also occurs near an N c = 0. 75 for a near optimum 
tube inside diameter of 5/8 inch. 

Plotting the maxima of figures 10(a) and (b) on figures 11(a) and (b), re- 
spectively, and the results for additional inside tube diameters yield the 
maximum value of Qrej/^ and its corresponding inside tube diameter for the 

three values of 5 s /& a chosen for the comparison. The curves given in figure 
ll(a) for the 500-kilowatt case indicate that maximum Qrej/^ for any choice 
of 5 s /& a occurs at a tube inside diameter between 3/8 and l/2 inch with the 
Qrej/ W curve relatively flat in this region. Maximum Q re j/W occurred at a 
diameter between l/2 and 5/8 inch for the 1-megawatt system as shown in fig- 
ure ll(b). 

The weight results of the double fin-tube configurations are compared to 
the central fin-tube results in figures ll(a) and (b) in order to show the de- 
sirability of the double fin-tube with reduced tube side-wall thickness. It is 
observed from figure 11 that the central fin-tube configuration has a larger 
heat rejection per unit weight than the double fin- tube with 5 s /5 a = 1.0. For 
the 500-kilowatt case, the double fin tube with 5 s /6 a =0.5 affords an 
11-percent weight advantage at maximum Qrej/w over the central fin tube, and 
a maximum upper limit of 39 percent when 5 s /5 a = 0. The radiator weight per 
kilowatt of electrical power for the 500-kilowatt case using a columbium radi- 
ator is 10.2, 15.0, and 19.0 pounds per kilowatt for the S s /S a = 0, 0.5, and 
1.0 cases, respectively. 

The percentage weight savings are reduced for the 1-megawatt power level 
with the beryllium radiator and are 8 percent for the 5 s /5 a = 0.5 case and 
32 percent for the 5 s /6 a = 0 case. For this power level, the specific weight 
was 3.50, 4.68, and 5.73 pounds per kilowatt for the 5 s /5 a = 0, 0.5, and 1.0 
cases, respectively. The beryllium radiator has a smaller percent increase in 
Qrej/W than the columbium radiator because there is less relative weight in 
the beryllium armor. 

The conclusions reached in figure 19 of reference 7, which compared the 
double fin-tube and central fin-tube conf igurations without taking into account 
the effects of header heat rejection and weight as well as pressure drops in 
the radiator tubes, showed a similar percentage increase in Q-rej/^ by reduc- 
ing & s /& a . This data, which was for a 1-megawatt, 1700° R, beryllium radia- 
tor, indicated a 6-percent weight advantage for the double fin tube at a 
5 S /S a =0.5 and a 28 percent advantage at 6 s /5 a = 0 over the central fin- 
tube configuration. 

The conductance parameter obtained at maximum heat rejection per unit 
weight is plotted as a function of inside tube diameter for a 5 s /5 a =0.5 Tor 
both geometries in figure 12. Results obtained for 5 s /5 a other than 0.5 
showed very little variation in the optimum value of N c . Both the 500- 
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kilowatt and 1 -megawatt examples are shown with the double fin-tube results 
yielding the higher value of H c for a specific choice of inside tube diam- 
eter. It is also noted that the values of H c obtained for the 500-kilowatt 
columbium radiator are somewhat larger than those obtained for the 1-megawatt 
beryllium radiator. Reference 7 indicated the double fin-tube configuration 
reached maximum Qrej/^ muc ^ larger values of N c than those obtained in 
this report. 


Radiator Geometry 

Flanform area . - The planform area of the double fin-tube geometry is 
shown plotted in figure 13 for the two power levels chosen for this investiga- 
tion. This figure illustrates the calculated variations of planform area with 
L/R 0 ratio, conductance parameter N c , and tube side -wall thickness to armor 
thickness ratio 6 s /& a for the two power cycles for peak Q re j/^ conditions. 
The two inside tube diameters used for the results shown in figure 13 corre- 
spond to near minimum weight conditions. Reducing the tube side-wall thickness 
ratio for peak L/r q ratios results in only a small variation in planform area 
for a constant value of N c but results in a sizable decrease in L/R 0 . How- 
ever, it should be kept in mind that the actual fin length Z is not decreas- 
ing that rapidly since the fin length is being increased by the amount 6 a - 5 S 
as 5 s /6 a approaches zero. The decrease in L/R q (which is a measure of the 
tube center-to-center distance) as S s /6 a decreases causes a reduction in tube 
block effectiveness (due to reduced Q^) and an increase in fin effectiveness 
(fig. 6) that result in only a small variation in total fin-tube effectiveness 
for a specific choice of N c . This accounts for the negligible variation of 
planform area with varying &s/^a- Similar results are obtained for the other 
inside tube diameters investigated. 

Comparison of the planform area results of the double fin-tube geometry 
with those of the central fin-tube geometry indicates good agreement as shown 
by figure 14. The planform area obtained for the central fin tube is less than 
that of the double fin tube for internal tube diameters greater than l/2 inch 
for the 500-kilowatt columbium radiator. At the 3/8-inch diameter, which cor- 
responds to near minimum weight, the two fin-tube geometries agreed to within 
4 percent. For the 1-megawatt beryllium system, the double fin tube has the 
larger planform area at tube diameters greater than 5/8 inch. At near optimum 
weight corresponding to a 5/8-inch diameter, the central fin-tube geometry 
offers only a small savings in planform area. The planform area obtained for 
the double fin- tube geometry in reference 7 for the 1-megawatt radiator is much 
greater than that obtained from the comparison calculations of this report. 

This is brought about because the conductance parameter at maximum Q re j/W for 

the case of reference 7 was larger than the values obtained in this report. 
Figure 13 verifies that increasing conductance parameter N c results in in- 
creased planform area. 

Fin thickness . - The effect of reducing the tube side-wall ratio 8 s /5 a 
on fin thickness is shown for the two examples in figure 15 for peak Q re j/W 
conditions at near optimum inside tube diameter. In both cases, reducing 
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& s /& a results in a sizable decrease in fin thickness at constant N c . Accom- 
panying this decrease in fin thickness is a reduction in L/R 0 and thus tube 
spacing. The magnitude of the fin thickness obtained for the two sample cases 
at optimum values of N c (fin thickness t is greater than 0.020 in. for the 
columbium radiator and greater than 0.035 in. for the beryllium radiator) re- 
sulted in values of fin thickness that are reasonable and should satisfy struc- 
tural and fabricational requirements. 

The comparison of the fin thickness obtained for the double fin-tube con- 
figuration indicates an increasing total fin thickness with increasing values 
of tube side wall to armor ratio 6 s /6 a throughout the range of Dj_ investi- 
gated for both power levels (fig. 16). Curves are given for fin thickness at 
maximum Q^j/W for both geometries and both system power levels with central 
and double fin-tube 5 s /S a =0.5 cases agreeing closely for the 1 megawatt 
radiator . 

Panel aspect ratio . - Radiator panel aspect ratio , which is defined as the 

ratio of panel width W to the tube length Z for the four panel radiator of 
figure 1, is shown plotted in figure 17 for the two double fin-tube example 
cases. The aspect ratio of a panel is seen to decrease as the tube block side- 
wall ratio & s /b decreases at a constant N c - The decrease in aspect ratio 
with reduced & s /6 a is primarily a result of the reduced l/R q since reduc- 
tions in & s /£> a had little effect on planform area as shown by figure 13. 

These results can be indicated by inspection of equation (48), which describes 
planform area and length, and equation (23), which describes the total panel 
width W. The reduction in aspect ratio at constant planform area is explained 
by the calculations which showed that as 5 s /s a gets smaller the number of 
tubes and the tube length increase as the L/R 0 ratio decreases . This was 
accompanied by a large reduction in the working fluid inlet tube velocity. 

The panel aspect ratio showed practically no difference between the cen- 
tral fin-tube and double fin-tube configurations as indicated in figures 18(a) 
and (b) for the 500-kilowatt and 1-megawatt systems, respectively. 

Number of tubes . - The number of radiator tubes required was found to in- 
crease substantially as the tube block side-wall ratio decreased. The 500- 
kilowatt columbium radiator with a tube inside diameter of 3/8 inch had 298, 
337, and 430 tubes for 5 s /S a ratios of 1.0, 0.5, and 0, respectively, at 
maximum heat rejection per unit weight. The 1-megawatt beryllium radiator 
using a tube inside diameter of 5/8 inch had 223, 250, and 302 tubes for 5 s /b a 
ratios of 1.0, 0.5, and 0, respectively, at maximum Q^j/W. 

For simplicity and reliability of fabrication it is desirable to reduce 
the number of tubes. A sizable reduction in the number of tubes without much 
reduction in Q^j/W can be had by just increasing the size of the tube inside 
diameter. For the 500-kilowatt columbium radiator with tube block side -wall 
ratio set at 0.5, a 3/8-inch inside tube diameter required 337 tubes whereas a 
l/2 -inch inside diameter results in 206 tubes. This reduction in the number of 
tubes results in only a 2-percent decrease in Q^, e j/w. l-^egawatt beryl- 
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lium radiator with 5 g /& a set at 0.5 required 250 tubes for a 5/8-inch inside 
tube diameter whereas a 3/4-inch inside diameter required just 182 tubes. This 
reduction resulted in only a 3-percent decrease in Qj-gj/W. Similar reductions 
in the number of tubes can be obtained at values of 5 s / 5 a equal to 0 and 1.0 
with small reductions in Q^gj/ W. 

Comparison of the number of tubes for the central and double fin-tube 
radiators at maximum Q^./w for the 500-kilowatt case indicated the central 
fin-tube geometry required 308 tubes compared to 298 tubes for the double fin- 
tube & s /S a =1-0 radiator at an inside tube diameter of 3/8 inch. This trend 
also held for the 1 -megawatt case for which the central fin-tube required 226 
tubes compared to 223 for the double fin-tube 8 s /& a =1.0 radiator. 

Interpolation of results. - If the values of the ratio of 60/&0 other 
than the three values chosen for this investigation (0, 0.5, and l) are desig- 
nated by forthcoming meteoroid bumper results, additional points may be ob- 
tained from curves faired through the calculated points to obtain intermediate 
values of panel aspect ratio W/Z^ fin thickness t, radiator planform area 
Ap, and radiator heat rejection per unit weight Q re j/tf. 


SUMMARY OF RESULTS 

An analysis of the double fin-tube configuration with variable tube block 
side -wall thickness and a comparison with a comparable central fin-tube radia- 
tor for the sample cases of a 1 -megawatt beryllium radiator and a 500-kilowatt 
columbium radiator for a 500-day mission time and probability of no puncture 
P(0) of 0.995 showed that 

1. A substantial increase in radiator heat rejection per unit weight is 
possible for the double fin-tube geometry if the tube block side-wall thickness 
can be reduced to less than 75 percent of that required for normal armor thick- 
ness . 


2. The double fin-tube configuration with ratio of tube block side-wall 
thickness to normal armor thickness of 0.5 offers an 11 and 8 percent weight 
savings over the central fin-tube geometry for the 500-kilowatt and 1-megawatt 
cases, respectively. This corresponds to 15.0 pounds per kilowatt for the 500- 
kilowatt case and 1.68 pounds per kilowatt for the 1 -megawatt case considered 
herein compared to 16.7 and 4.83 pounds per kilowatt for the central fin-tube 
500-kilowatt and 1 -megawatt cases, respectively. 

3. The percent weight savings is increased to 39 and 32 percent for the 
500-kilowatt and 1 -megawatt cases, respectively, if the ratio of tube block 
side -wall thickness to normal armor thickness is set at 0. In this case, the 
corresponding radiator specific weights would be 10.2 and 3.50 pounds per kilo- 
watt for the double fin-tube geometry and 16.7 and 4.83 pounds per kilowatt for 
the central fin-tube geometry. 

4. The double fin- tube configuration reached maximum heat rejection per 
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unit weight for the 1-megawatt system at a conductance parameter of 0. 75 where- 
as maximum heat rejection per unit weight for the central fin-tube geometry 
occurred at 0.63. These values for the 500-kilowatt case were 0.63 and 0.75, 
respectively, for the double and central fin-tube. However, the conductance 
parameter associated with maximum heat rejected per unit weight decreases with 
decreasing inside tube diameter for both the double and central fin-tube radi- 
ators. 


5. The double fin-tube radiator reaches maximum heat rejected per unit 
weight for a tube inside diameter in the range from 3/8 to l/2 inch for the 
500-kilowatt case, and from l/2 to 5/8 inch for the 1-megawatt case. 

6. The physical dimensions of the double fin-tube radiator (planform area, 
panel aspect ratio, number of tubes, inside tube diameter, and fin thickness) 
can be varied over a fairly wide range without seriously decreasing the radia- 
tor heat rejection per unit weight. 

7. The fin thickness obtained for maximum heat rejection per unit weight 
for the double fin-tube geometry is of reasonable fabricational and structural 
magnitude (greater than 0.02 in. for the 500 kw columbium radiator and greater 
than 0. 035 in. for the 1 Mw beryllium radiator) . 

8. The values of planform area and panel aspect ratio obtained for the 
double and central fin-tube geometries were in close agreement regardless of 
the choice of tube block side-wall ratio near the minimum weight condition. 
Comparison of the total fin thickness for the two geometries agreed well at a 
tube block side -wall thickness ratio of 0.5, but showed considerable varia- 
tions at side-wall ratios of 1.0, and 0. 

9. Reduction of the tube block side- wall thickness ratio substantially 
increases the number of radiator tubes at maximum heat rejection per unit 
weight. However, the number of tubes can be reduced to the values correspond- 
ing to the central fin-tube geometry by increasing the value of Dj_ above 
that for maximum heat rejection per unit weight. The weight penalty involved 
is only several percent. 


Lewis Research Center 

Rational Aeronautics and Space Administration 
Cleveland, Ohio, August 12, 1964 
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Figure 3. - Central fin- tube geometry. 
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Figure 4. - Assumed radiator heat- transfer paths. 
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Figure 7. - Tube effectiveness; double fin-tube 
geometry. Powerplant output, 500 kilowatts; 
temperature at radiator inlet, 1700° R; in- 
side diameter, 3/4 inch; columbium radiator. 
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(b) Powerplant output, 1 megawatt; inside tube diameter, 5/8 inch; beryllium 
radiator. 

Figure 9. - Heat rejection per unit weight characteristics; double fin-tube 
geometry. Temperature at radiator inlet, 1700° R; ratio of tube side-wall 
to armor thickness, 0.5. 
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Maximum heat rejection per unit weight, Q/w, Btu/(hr)(lb) 
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(a) Powerplant output, 500 
kilowatts; columbium 
radiator. 


(b) Powerplant output, 1 megawatt; beryl 
Hum radiator. 


Figure 11. - Maximum heat rejection per unit weight. Temperature at radiator 
inlet. 1700° R. 
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Figure 12. - Conductance parameter at maximum heat rejec 
tion per unit weight. 
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(a) Powerplant output , 500 kilowatts; columhium radiator. 
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(b) Powerplant output, 1 megawatt; beryllium radiator. 


Figure 14. - Planform area comparison at maximum heat re- 
jection per unit weight for double and central fin-tube 
radiators. Temperature at radiator inlet, 1700° R. 
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(a) Powerplant output, 500 kilowatts. 




(b) Powerplant output, 1 megawatt. 

Figure 18. - Panel aspect ratio comparison at maxi- 
mum heat rejection per unit weight for double and 
central fin-tube radiators. Temperature at 
radiator inlet, 1700° R. 
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